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Abstract 

We discuss high energy hadronic collisions within the theory of the Color Glass Condensate. 
We point out that the initial electric and magnetic fields produced in such collisions are longi- 
tudinal. This leads to a novel string like description of the collisions, and a large Chern-Simons 
charge density made immediately after the collision. The presence of the longitudinal magnetic 
field suggests that essential to the description of these collisions is the decay of Chern-Simons 
charge. 
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1 Introduction 



High energy hadronic collisions can be described as collisions of sheets of Colored Glass 
Condensate [1,2,3,4,5,6,7,8,9]. The degrees of freedom of the Color Glass Condensate are 
those of high energy density gluonic fields [10,11,12,13,14,15,16]. Because the physical den- 
sity of gluons becomes large, the typical separation between gluons is small, and therefore 
a s is small. The highly coherent gluons fill the phase space up to the maximal occupation 
number ~ l/a s and can thus be thought of as condensed. Because of their high speed 
and Lorentz time dilation, the valence degrees of freedom are seen by the low x fields as 
slowly evolving in lightcone time. Systems which evolve over long time scales compared 
to natural ones are glasses. Hence the name Color Glass Condensate. It is possible to 
describe this system from first principles in QCD because the coupling is weak. 

The equations which describe the color fields in a high energy hadronic collision have 
been written down in Refs. [1,2,3,4]. For the case where the rapidity dependence of 
the fields can be ignored, the equations were numerically solved [5,6,7,8,9]. A success- 
ful phenomenology has been applied to RHIC energy collisions of nuclei and dA col- 
lisions [17,18,19,20,21,22,23,24,25,26,27]. Since then, a literature has developed relating 
this description to issues such as topological charge generation [28,29,30,31] and plasma 
instabilities [32,33,34,35,36,37] . 

Topological charge generation is related to a number of non-perturbative phenomena in 
field theory. In electroweak theory, it leads to the anomalous generation of baryon plus 
lepton number. In QCD, it is related to the violation of U(l) chiral symmetry, which 
may ultimately drive the breakdown of the SU(2) chiral symmetry and be responsible for 
mass generation in QCD. As the topological charge is CP odd, it may also make large 
CP violating fluctuations in heavy ion collisions [28,29,30,31]. Some indications of these 
fluctuations have been recently observed by the STAR expertiment [38] . 

The instabilities may generate the rapid thermalization seen at RHIC. This intermediate 
matter is highly coherent, and makes the transition from the Color Glass Condensate 
to the Quark Gluon Plasma. We shall call it the Glasma. This is because the fields 
are coherent and have strength 1/g, and so in their interaction with particle degrees of 
freedom, are 0(1), since the g of interacting with the field is canceled by the 1/g from the 
strength of the field. This enhances the rate at which the system can rearrange in phase 
space and reach a thermalized distribution. 

In this paper, we will discuss some properties of the Glasma. Our intent is not to add 
anything to the existing mathematical literature concerning the solution to the classical 
field equations resulting from the collision. In fact, we largely restrict ourselves to the 
somewhat tame situation where the fields are boost invariant, and we only briefly discuss 
plasma instabilities. 

Our goal is to provide a qualitative and intuitive understanding of the Glasma fields 
shortly after collisions. This work is an elaboration and interpretation of some of the 
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recent observations of Fries et. al. [39,40]. The structure we illuminate is amusing. An 
infinitesimal time after the collision, longitudinal electric and magnetic fields are produced. 
There are no transverse fields, except on the valence sheets of charge which are the sources 
for the Glasma field. In fact, infinitesimally after the collision, the charge distribution on 
these sheets is modified by the addition of sources of color electric and color magnetic 
charge. (We will discuss more precisely what is the physical approximation which leads 
to this infinitesimal change, and what is realistic physical time interval involved.) 

The physical situation immediately after the collision bears close analogy to string models 
of high energy collisions [41,42,43,44]. In models such as that advocated by the Lund 
group, there is a longitudinal electric but no longitudinal magnetic field. It is the decay of 
these fields which is the essential dynamics of the Glasma, and which ultimately produces 
the Quark Gluon Plasma. For the Glasma, these fields can decay due to both classical 
rearrangement of the field, or from quantum pair creation. The classical rearrangement of 
the field into radiation of gluons with pt ~ Q s is a highly coherent process and happens 
on the time scales typical of the inverse separation of the color charges. The quantum 
process should be a factor of a s weaker and naively appears to a be a small correction 
to the classical process. We shall argue that these processes are related. We will discuss 
the recent proposal by Kharzeev, Levin and Tuchin that such fields are precisely the type 
which yield exponential distributions in transverse mass [45,46]. 

Finally, a longitudinal electric and magnetic field has a non-zero topological charge den- 
sity FF. We show how this results in a non-zero Chern-Simons charge density. Such a 
density can seed interesting non-perturbative phenomena associated with chiral symmetry 
breaking. We also argue that the decay of these fields bears a remarkable resemblance to 
the picture advocated by Kharzeev, Kovchegov and Levin [30] and by Janik, Shuryak and 
Zahed [47,48], which interprets high energy nuclear collisions as the decay of instantons 
(see also Refs. [49,50,51] for a discussion of instantons in deep inelastic scattering). 



2 The Classical Equations 

For simplicity, we consider collisions at very high energy where the rapidity dependence of 
the distribution of produced particles can be ignored. We assume the density of produced 
particles is so large that the typical transverse separation of the produced gluons is very 
small compared to a Fermi. The typical QCD interaction strength is therefore a s < 1. 

We separate the degrees of freedom of a fast moving nucleus into large x and small x 
degrees of freedom. This separation is arbitrary, but leads to an effective Lagrangian for 
the small x degrees of freedom. The parameters of this Lagrangian are subject to renor- 
malization group evolution which ultimately resolves the ambiguity of the arbitrary scale 
separation, and determines the parameters of the effective Lagrangian. The fast degrees 
of freedom are Lorentz contracted and on a sheet, and produce fields which ultimately 
describe the low x degrees of freedom. 
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The fields of the single nucleus are the Lienard-Wiechert potentials associated with the 
color charge distribution. They exist also only in the sheet and for each source of charge 
are mutually orthogonal electric and magnetic fields which are also perpendicular to the 
beam direction. It would appear that there is a paradox since the quanta associated with 
these fields are at small x, although the spatial extent of the fields is in the region of the 
valence quanta with a longitudinal size which is very small. By the uncertainty principle, 
it would seem that the gluons associated with this field are at large x, not small. This 
paradox is resolved by the fact that the small x quanta are those of the Fourier transform 
of the vector potential (in lightcone gauge). Although the color electric and color magnetic 
field exist only within the sheet, the vector potential has a discontinuity at the sheet but 
exists in the region outside the sheet. This is where the small x gluons are located, so 
that the large longitudinal spatial extent indeed corresponds to small x. 




Fig. 1. Color fields in spacetime. In regions (1) and (2), where only one of the nuclei has passed 
by, the field is the pure gauge field of this one nucleus. In region (3) the field is known numerically. 

The sheet of color charge is approximated as infinitesimal, but its physical size can be 
estimated by renormalization group arguments. Let us concentrate on the nucleus moving 
along the x~ = Odightcone. There is a cutoff A + separating the hard and soft degrees of 
freedom. The degrees of freedom with p + > A + are described as classical color charges 
and the ones with k + < A + as classical fields, characterized by the saturation scale Q s . 
The renormalization group equation for this system describes how the distribution of 
color charges evolves when the cutoff is changed. When A + changes by a factor e Ay , the 
correction to the saturation scale is of order a s Ay and becomes of order one only when 
Ay ~ l/a s . When studying the soft degrees of freedom around midrapidity y ~ r\ ~ 
we can therefore assume that they are separated from the sources with p + > A + by a 
rapidity difference Ay ~ l/a s . The characteristic momentum and time scales for these 
soft modes at y ~ are k + ~ Q s and At ~ 1/Q S . These soft classical fields see hard 
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sources as localized in an interval Az ~ Ax~ ~ l/p + < z~ 1 ^ as /Q s which, being close to 
the x~ = O-lightcone (see Fig. 1), corresponds to a time interval At ~ e~ 1 ^ (Xs /Q s . The 
smallness, in the weak coupling limit, of this scale e~ 1/,as /Q s compared to the scale 1/Q S 
of the soft fields is the physical picture approximated by the infinitesimally thin sheet. 

For the collision of two nuclei, we imagine working in the center of mass frame. Here 
the two nuclei can be thought of as sheets of colored glass. The fields associated with the 
collision are illustrated in Fig.l. We choose the fields to vanish in the backward light cone. 
On the side light cones, we choose the fields to be two dimensional pure gauge transforms 
of vacuum, 

4,2) = Y g U ^ U l2) (2- 1 ) 

The derivatives are two dimensional and transverse in the sheets. The discontinuity of the 
field at x ± = in the backward light cone is determined by the source density on the two 
sheets, 

r = 5<* + 5(x-) Pl (yL T ) + 5<*-5(x + )p 2 (x T ) (2.2) 

In practice, to relate the vector potential to the charge density, one has to spread the 
charge density out in the longitudinal space, but this does not affect the determination of 
the solution of the classical field in terms of the vector potentials A l ^ 2 y 

In Fig. 2, the classical color and magnetic fields are shown on the two nuclei before the 
collision. They are frozen in time and the color electric and color magnetic fields are 
orthogonal to the beam direction and to one another. 

In the forward light cone, there is no solution which is pure gauge. There is therefore 
particle production and evolution in the forward lightcone. This matter is the Glasma, 
until it eventually evolves into a thermalized Quark Gluon Plasma. Because the sources, 
Eq. 2.2, are restricted to the light cone there is a boost invariant solution of the field 
equations in the forward light cone. It is convenient to work in the r, 77 coordinate system 
and in radial gauge, 

A T = i (x + A~ + x~A + ) = 0. (2.3) 

The field can then be written in terms of 1 a = A v and the transverse components as 
A ± = ±x ± o;(r, xt) 

A^afaxT). (2.4) 

The equations of motion for these fields were written explicitly in [1,2]. As second order 
partial differential equations they require initial conditions for the fields and the first time 



1 The Hamiltonian equations of motion for the numerical calculation are often, e.g. [5,9] written 



in terms of (f> = A v = — r 2 ^ 7 . 
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Fig. 2. The fields on the two nuclei prior to the collision. 

derivatives. These boundary conditions at r = are determined by requiring that the 
Yang Mills equations be solved across the forward lightcone, and are 



a\ | T=0 = a\ + a\ 



a |r=o 
d T a | T=0 



2 

d T a 



3 k=o= 



0. 



(2.5) 



If one then solves the Yang-Mills equation near the light cone, one finds that the transverse 
color electric and color magnetic fields vanish as r — > 0, but that the longitudinal electric 
and magnetic fields 2 are non- vanishing [39,40], 

E z = ig[a{,ai] 

B z = tge^[a\,ai]. (2.6) 

(2.7) 

A plot of the transverse and longitudinal color fields is shown in Fig. 3, based on a numerical 
solution of these equations. 

A picture of the fields after the collision is shown in Fig. 4, at a time infinitesimal after the 
2 In the sense of the usual t, ^-coordinates, i.e. E z = F tz and B z = F xy = \^F % K 
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Fig. 3. Components of the gauge field, computed numerically on a 512 2 -lattice with 
g 2 »R A = 67.7. 

collision. In addition to the fields in the sheet, there are longitudinal electric and magnetic 
fields, which are random in color with a correlation length scale in the transverse plane of 
order 0(1/ Q s ). There are new sources of color electric and magnetic charge on the sheets 
on which the original fields resided. 

The presence of the longitudinal fields is at first sight surprising. How did the longitudinal 
electric and magnetic field arise? Let us define the chromoelectric and chromomagnetic 
charge densities by 

V • E — Pelectric 

V • -B Pmagnotic- (2-8) 

Because the fields satisfy the nonabelian Gauss law and the Bianchi identity these charge 
densities have a contribution from the gauge field itself: 

pelectric =ig{A\E 1 } 

Pm ag netic=i9[A\B i ]. (2.9) 

The source terms turn on instantaneously in the forward light cone. This arises because 
there are transverse electric and magnetic fields in the sheets, which in the forward light 
cone are multiplied by a no n- vanishing vector potential from the opposing sheet. The 
induced charge density on one sheet is the negative of that on the other, as it must be 
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Fig. 4. The color electric and magnetic fields after the collision. In addition to the transverse 
fields on the sheets there is a longitudinal field between them, originating from electric and 
magnetic charges on the sheets. 

to generate the longitudinal magnetic and electric fields and be consistent with Gauss's 
law. To be explicit, the nucleus on the x~ = light cone has a gauge potential A 1 = 
9{x~)a\{xT) and the other nucleus A 1 = 9{x Jr )ot\(x.T)- The cross terms of these two give 
rise to the opposite charges on the light cones, electric 



Pelectric 



V2 



9(x )5(x + )[a\, a\ 



W_ 
V2 



9(x + )5(x~)[ai,a[] 



and magnetic 



(2.10) 



Pmagnctic = ^ 6 ( X ~ ) 5 ( X+ ) ^ H > <4] + ^f^ 6 ^ ^ M ' ^1 , 

which exactly correspond to the initial longitudinal fields, Eq. (2.6). 



(2.11) 



After the collisions, sources of magnetic and electric charge have appeared on the sheets. 
The valence structure of the nuclei have changed, and there are new sources for the 
fields. This description is very similar in spirit to the original suggestions for a string 
based description of heavy ion collisions. These old pictures had only sources of electric 
charge [41,42,43,44]. The appearance of a longitudinal magnetic field is an entirely new 
aspect. 

In the Lund Monte Carlo model [41], there are longitudinal electric fields induced by a 
collision, which subsequently decay by quantum pair production. The situation here is 
similar except for the longitudinal magnetic field and the different length scale in the 
transverse direction, 1/Q S in stead of 1/Aqcd- The other essential difference is that the 
fields can decay away by classical evolution of the charged Yang-Mills field. There is not a 
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necessity for pair production from quantum effects to quench the fields. We will expand on 
this later, since, the quantum and classical quenching processes are related in a non-trivial 
way. 

It is interesting to note the structure of the energy momentum tensor T^ v = \g^ u F af5 F a p — 
F^Fva for this initial condition. It is diagonal and, as always in gauge theory, traceless: 
Tftv = \{F>\ + Bl) x diag(l, —1, — 1, 1). This can be compared to the standard form for 
a system with an anisotropy in the ^-direction: T^ v = diag(e, — p±, —p±, — Pl), where e is 
the energy density and p± and pi are the transverse and longitudinal pressures. We see 
that the initial field configuration has negative "longitudinal pressure" . The configuration 
that is the starting point for studies of isotropization by plasma instabilities, where the 
diagonal elements of T^ v at rj = are (e, — e/2, — e/2, 0), is only reached at times r > 1/Q S 
when the classical fields start to behave linearly due to the expansion of the system. 

Finally, the non-zero longitudinal electric and magnetic fields implies that the initial 
conditions have a large density of topological charge F^F^. Since 

d^K, = F^F^ (2.12) 

there is also a nonzero Chern-Simons current. We will investigate the topology of these 
configurations in the next section. The physical picture we have of the initial conditions 
looks quite similar to that proposed by Kharzeev, Kovchegov and Levin [30] and by Janik, 
Shuryak and Zahed [47,48], who claimed that high energy collisions might be described by 
the decay of instanton like configurations, the sphalerons. The decay of these longitudinal 
electric and magnetic fields is in fact the decay of topological charge. 

To thermalize the system, one must first produce quanta associated with the gluon field. 
In order to be described as quanta, the classical field associated with these quanta must 
be weak, A <C 1/g. For a mode with momentum p?, this occurs in a time of order 
t ~ l/pr- For a time l/pr <C r <C l/(a s pT), the field has strength 1 « A « 1/g, 
and can be treated either classically, or as an ensemble of particles. This follows because, 
up to oscillations, the field behaves as 1/ ' y/r at late times. For p T ^> l/(asr), the field 
is essentially quantum, as the classical field is small compared to the field induced by 
quantum fluctuations (see [52] for an argument approaching this same limit from the 
opposite direction, i.e. kinetic theory, and [53,54] for a more formal argument). 

Roughly speaking, the emission of particles occurs during the time interval 1/Q S <C r <C 
l/(a s Q s ). Since a s varies logarithmically with Q s , the time scale is crudely l/Q s - Nev- 
ertheless, this treatment suggests there is an intermediate time scale where either the 
classical and transport solutions to the evolution of a coupled system of fields and par- 
ticles may be a good approximation. Above the high momentum end of this matching 
interval, one should use a hard particle treatment. Below the lower scale of the overlap 
region, one should use a classical field treatment. For a complete treatment, one should 
have a floating scale which depends upon time. Modes above this floating scale are hard 
particles, and those below are soft. 
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The picture is therefore of classical fields evaporating into gluons, and a system where the 
hard gluons are interacting with a classical field. The interaction of the hard gluons among 
themselves is probably not so important during the time interval 1/Q S <C r <C l/(ot s Q s ) 
since the characteristic time scale for hard particles to thermalize by interactions among 
themselves is of order r ~ l/(a^Q s ), since these interactions involve scattering cross 
sections. The interaction of the classical field with itself is important, and the characteristic 
time scale for interaction of the hard fields with the classical fields is also important, since 
the effects of coupling constant cancel out in the interactions of the hard field with the 
coherent fields strength A ~ 1/g. 

It is this system which we refer to as the Glasma. 

An issue of current interest is whether or not the boost invariant solution described in 
the previous section is stable with respect to space-time rapidity dependent fluctuations. 
There is reasonably strong evidence that the coupled particle-field Glasma system has 
such instabilities. The issue of the purely classical evolution also has claims that this may 
be the case [37], but they happen over extremely long time scale, and one may be worried 
about whether such instabilities reflect the continuum limit. 

As we will discuss in Sec. 3, the topological properties of the gauge fields depend cru- 
cially on the dimensionality of the system. One would therefore expect that relaxing the 
strict boost invariance of the field configurations will significantly change the dynamics of 
the Chern-Simons charge. Perhaps the instability of the boost invariant Glasma and the 
production of net topological charge might be related. 

If the purely classical system is unstable, then one may legitimately worry about the 
validity of the entire Glasma description. The problem is that if the instability is like 
that of classical turbulence, then different initial conditions which are separated from one 
another infinitesimally, will at late times evolve to solutions far apart with AS* ~ S . Such 
a difference in initial conditions might be generated by quantum noise. Their contribution 
to the path integral will differ by a factor of 

e iAS (2.13) 

which is oscillatory. Therefore one would not expect convergence of the glassy sum over 
initial conditions. 

There can also be instabilities with AS ~ a s S <C S. In addition to the interpretation 
of the momentum modes of the classical field as gluonic quanta there is quantum radi- 
ation. This may be generated by considering the small fluctuations in the classical time 
dependent background field. The effective action is then obtained by integrating over 
these fluctuations. At fourth order in the expansion of the effective action in powers of 
the external field, there is an imaginary part generated, coming from diagrams such as 
the one in Fig. 5. This imaginary part is precisely pair production. Pair production does 
not correspond to an instability of the purely classical equations, but to an instability of 
the classical background field in a quantum mechanical treatment. It remains an effect 
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Fig. 5. Left: A lowest order (both in the classical field and in the quantum corrections) diagram 
that contributes to the imaginary part of the effective action. The winding lines represent the 
classical field, 0(l/g) and the wavy lines the quantum fluctuation, 0(1). Right: The correspond- 
ing pair production diagram. 

suppressed by a s with respect to the classical background and is therefore not a problem 
for the validity of the glassy description. A consistent computation of this pair production 
must also include the renormalization group evolution of the classical fields, which con- 
tributes at the same order in a s [55]. Kharzeev, Levin and Tuchin [45,46] have argued that 
for fields of the type we have shown exist in the collision of two sheets of colored glass, 
specifically a pulsed longitudinal electric field, the distribution of produced particles is an 
exponential in transverse mass, contributing to the early thermalization of the system. 

The contribution to the gluon multiplicity from pair production is related to the diagram 
shown in Fig. 5, where one integrates over one of the legs of the quantum fluctuation 
field. In this integration the integral over Ay, the rapidity gap between the produced 
gluons, diverges (see e.g. Appendix B of [56]). The region Ay > l/a s corresponds to 
renormalization group evolution of the sources and the genuine quantum correction to 
gluon production only comes from the region Ay <C l/a B . To compute the first quantum 
correction to gluon production one must introduce a cutoff Y ~ 1 /a s in the integration 
over Ay. One then integrates over one of the legs in Fig. 5 up to this cutoff. The physical 
result is independent of Y, because the dependence of the gluon multiplicity on Y must 
be cancelled by the renormalization group evolution of the sources as a function of Y. 
Thus pair production remains suppressed by a power of a s compared to the background 
field. This can be contrasted with the calculation of quark pair production [57,58,59,60], 
where the integral over Ay is convergent and the quark multiplicity is directly suppressed 
by a s compared to the classical field. Quark pairs, unlike gluons, do not contribute to the 
renormalization group evolution of the source to this order. 



3 Topological Charge 

The topological charge density of a gauge field configuration is defined as 
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d v a a — - 

— = -^TtF^F^ = -^TrE ■ B, (3.14) 
da; 4 32tt 2 8tt 2 k j 

where F^ v = \e^ V(>C! F pa . The charge density is the four- divergence of the topological 
current, F£ v F pv = d^K^, where 

K" = e^ pa A a u (V; CT + ^f abc A b p A c ^ . (3.15) 

Integrating the 0-component of the current over space one can define a topological charge 
v = g£p- / d 3 xfT°. The change in the topological charge between times, 

v{t -> oo) - u(t -> -oo) = I dtdtv = I d 4 zTrF^„, (3.16) 

is a gauge invariant quantity, but the topological current and charge themselves are gauge 
dependent. It is therefore essential in our case that one cannot transform away the pure 
gauge contribution at large r, because such a gauge transformation would introduce in 
the backward lightcone a nonzero gauge potential A^. 

Working in the boost invariant case it is more natural to look at the Chern-Simons charge 
per unit rapidity on a constant proper time surface 

% = ^l d2xTTKr = 64^ / d2xre " K (*3 + 9f abc AtA<) + 2A°F« V ] (3.17) 

As was pointed out in [31] this expression can be written in a much simpler form when the 
field configurations are boost invariant. Performing a partial integration one can rearrance 
Eq. (3.17) as 

The surface term from the partial integration is 

-.2 ,. Jl 
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327T 2 



J d 2 x T e^- (A* A;) = J— 2 j>MA'}A a r (3.19) 



There should be a large topological charge density associated with the Glasma. This is 
because the Glasma fields at r — > + have longitudinal E and B and give nonvanishing 
E • B. Of course, the typical correlation length in these fields in transverse size is the sat- 
uration momentum, and for reasons described below, we expect zero total Chern-Simons 
charge from these fields. It is interesting that this resembles the situation envisioned by 
Kharzeev, Kovchegov and Levin [30] and Janik, Shuryak and Zahed [47,48] in that the 
production of particles can be thought of as due to the decay of Chern-Simons charge. 
Presumably, this is somewhat like the situation at finite temperature where one has a 
density of sphalerons associated with the decay of topological charge. It is amusing here 
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that the topological charge density arises already at the classical level and is 0(1), not 
0(a a ) as it would be from quantum fluctuations. 

For early times the electric and magnetic fields are both longitudinal, and there is a 
nonzero topological charge density. Note that because d/di] = rd/dz the charge density 
per unit rapidity is zero, but the density per unit z nonzero. The expectation value of 
the Chern-Simons charge density when averaged over the source color charge densities is 

zero, but the magnitude of the fluctuations, (^(du / dr/) 2 ^ , can be estimated as follows. The 

system consists of approximately N = nR 2 A Q 2 uncorrelated domains, each having charge 
~ 1. The sum of these independent charges will fluctuate around zero with magnitude 



(du/dr]) J ~ yN ~ Q s Ra- This estimate is confirmed by the numerical calculation 
[31]. 

For late times r — > oo the field gets weak due to the expansion of the system and the gauge 
field can be written as a gauge transformation of a solution to the linearized equations of 
motion [2]. The general such asymptotic form of the solution is: 

Ai = V(x. T ) (q(x t , r) - l -d)j V\x T ) 

A, = \/(x T )e(x T ,r)\/ t (x T ) (3.20) 



with 

ei(*r,T)=f l^ 1 ™ [j (\k T \r)h{(k T ) + N (k T r)hf(k T ) 

e(x T , r) =J l^ 1 ™ r [j^r)^ (k T ) + iV 1 (k T r)/ i JV (k T )] . (3.21) 

This decomposition between the field and the pure gauge contribution is not unique, but 
this ambiguity can be removed by requiring that q satisfies the transversal Coulomb gauge 
condition <9j6j = 0. The mode functions hj (k T ), hf(k T ), h J (k T ), h N (k T ) are determined 
by the nonlinear dynamics of the system at early times. 

In the lowest order perturbative solution [2] the mode functions are determined only by the 
initial conditions. Because the Neumann functions N n diverge at the origin, hfikr) and 
h N (k T ) vanish at in the lowest order. Thus the Chern-Simons charge, Eq. 3.18, for each 
transverse momentum mode behaves as J (|k T |r) J^k^r), giving a zero contribution 
when averaging over a time period r ^> l/|k^|. At higher orders in the source also the 
Neumann function solutions can be present and hf ikx) and ^(k^) can be finite. Because 
Jo(|kr|r)A r i(|kr|r) and iV (|kr|r) Ji(|k T |r) oscillate around finite values, these terms can 
give a finite radiative contribution to the topological charge. This higher order radiative 
contribution is the most likely interpretation for the finite topological charge seen in the 
numerical calculation [31]. 

When, as in Ref. [31], the topological charge is computed on a finite lattice with periodic 
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boundary conditions and thus the surface term Eq. (3.19) vanishes exactly. In general it is 
possible to have a solution where the combination of a covariantly constant longitudinal 
field A v = V^x^eV^xy) with d r e = die = and the pure gauge component V(x. T ) give 
a nonvanishing topological charge. This contribution is not present in the perturbative 
solution of Ref . [2] , because a constant A v corresponds to A v that diverges at r = and 
is not allowed by the initial condition. For an arbitrary color charge distribution in the 
transverse plane the pure gauge fields of the individual nuclei, and thus also the pure 
gauge fields in the asymptotical solution Eq. (3.20), vanish only logarithmically for large 
x. T . However the overall color charge of the nucleus must be zero, and thus for distances 
larger than 1/Aqcd outside the nucleus the color field must die off faster, meaning that 
the boundary term will vanish. 



3.1 Homotopy Groups 

The nth homotopy group 7i n (G) of a topological space G, in our case the gauge or sym- 
metry group, is the group of mappings from the n-sphere S n to G. The usefulness of 
homotopy groups in field theory often arises in the following way. Let us presume were 
are interested in a field theory with symmetry group G. When one requires that field 
configurations in R n approach a constant at large distances, one can view space as the 
compact group S n . The homotopy group n n (G) then tells us whether topologically in- 
equivalent field configurations exist. One example is the group 7^(0(3)) = Z, which 
tells us that in an 0(3)-symmetric nonlinear sigma model in 2 dimensions one can con- 
struct a topological charge that takes integer values [61]. Another example is the Skyrme 
model [62,63], where the solitons of an SU(2) symmetric "pion" field in three dimensions 
can be interpreted as having integer baryon number, because 7r 3 (SU(2)) = Z. For this 
same reason the Chern-Simons number of SU(3) in a gauge theory changes by an integer 
in a nonperturbatively large gauge transformation. The second homotopy group of SU(3), 
however, is trivial: 7T2(SU(3)) = 0, meaning that for boost invariant two dimensional field 
configurations one can not construct a topological charge taking integer values. 



4 Conclusions 

We have in this paper tried to elucidate some known but little appreciated aspects of 
Glasma, the highly coherent matter making in transition from Color Glass Condensate 
to Quark Gluon Plasma. It is characterized, in a manner reminiscent of the Lund model, 
by the decay of a longitudinal electric and magnetic field into particles. This decay can 
happen both as classical radiation of the field and, quantum pair production from the 
classical background. We have discussed the relation of these two processes and pointed 
out that a full computation of the pair production must be done consistently with the 
renormalization group evolution of the sources. 
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As the magnetic and electric fields are initially parallel, they correspond to a nonzero 
Chern-Simons charge density. The evolution of the Glasma can thus be described also as 
the decay of this Chern-Simons charge. 
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